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Abstract 

We present the geometric solutions to some variational problems 
of statistical mechanics and combinatorics. Together with the Wulff 
construction, which predicts the shape of the crystals, we discuss the 
construction which exhibit the shape of a typical Young diagram and 
of a typical skyscraper. 
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1 Introduction 



In this review we will discuss the phenomenological geometrical construction, 
called Wulff construction. Its importance lies in the fact that its predictions 
turn out to be true theorems of mathematical physics and combinatorics. 
Below we present some of the rigorous results in which this construction 
yields the answer to important questions. 

*The work was partially supported by the Russian Fund for Fundamental research 
through grant 99-01-00284. 
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1.1 Statistical mechanics 



The variational problems of statistical mechanics we are going to discuss here 
are those related to the formation of a droplet or a crystal of one substance 
inside another. The question here is: what shape such a formation would 
take? The statement that such shape should be defined by the minimum 
of the overall surface energy subject to the volume constraint was known 
from the times immemorial. In the isotropic case, when the surface tension 
does not depend on the orientation of the surface, and so is just a positive 
number, the shape in question should be of course spherical (provided we 
neglect the gravitational effects). In a more general situation the shape in 
question is less symmetric. The corresponding variational problem is called 
the Wulff' problem. Wulff formulated it in his paper |^] of 1901, where he 
also presented a geometric solution to it, called the Wulff construction (see 
section 2.2 below). 

This Wulff construction was considered by the rigorous statistical mechan- 
ics as just a phenomenological statement, though the notion of the surface 
tension was among its central notions. The situation changed after the ap- 
pearance of the book ||DKS|| . There it was shown that in the setting of the 
canonical ensemble formalism, in the regime of the first order phase transi- 
tion, the (random) shape occupied by one of the phases has asymptotically (in 
the thermodynamic limit) a non-random shape, given precisely by the Wulff 
construction! In other words, a typical macroscopic random droplet looks 
very close to the Wulff shape. The results of the book ||DKS|| are restricted 
to the 2D Ising ferromagnet at low temperature, though the methods of the 
book are suitable for the rigorous treatment of much more general models. 

Another problem where Wulff construction is of importance is the theory 
of metastability. This is a theory which studies the states of the matter 
which "should not be there" , but which still can be observed, albeit for only 
a short time. One example is water cooled below the zero temperature. 
This supercool water can stay liquid, but not for a long time, and it then 
freezes abruptly. Such states are called metastable. They should be thought 
about as continuations of the equilibrium states through the point of phase 
transition, at which these states cease to exist as equilibrium states. It turns 
out that their lifetime is again determined by the quantities given by the 
Wulff construction. 
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1.2 Combinatorics. 



Combinatorics is another field where similar problems appear. Specifically 
we have in mind the problem of finding the asymptotic shape of a typical 
Young diagram and related problems. This question was first asked and 
answered by Vershik and his coauthors, see |[VKei| , [Vl| , [V^ , PVZ]| . These 



papers were treating the two-dimensional case, while the higher dimensional 



problems were solved only recently, see |[Ke|| . In combinatorics there is also a 
variational problem to be solved, which in fact is close to the Wulff problem. 

When compared with the situation in statistical mechanics, the combina- 
torial program and its development look very similar. The only difference is 
that the combinatorial counterpart of the Wulff construction was not known, 
probably because there was no heuristic period there. Below we describe 
such a construction, presented in |^^. It provides, like the Wulff one, the 
geometric solution to the corresponding variational problem under minimal 
restrictions on the initial data, and also proves the uniqueness of the solution. 

The next section contains the necessary geometric statements. Some are 
given with proofs. In Section 3 we discuss the problems from statistical 
mechanics, in which these geometric problems appear naturally. The last 
Section 4 deals with the corresponding combinatorial problems. 



2 Geometry 

In this section we consider various geometrical problems. We explain the 
Wulff minimization problem and the Wulff construction (^), which solves it. 
We then describe a related maximization problem. We finish by discussing 
flat facets of the Wulff shapes. 

2.1 The droplet energy functional 

Let S"^ C M'^^^ denote the unit sphere, and let the real function r on S'^ be 
given. We suppose that the function is continuous, positive: r (■) > const > 
0, and even: r (n) = r (— n) . Then for every hypersurface M'^ embedded in 
gg^j^ define the Wulff functional 

Wr{M^) = [ r(n,) ds,. (1) 
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Here x G M'^ is a point on the manifold M'^, the vector is the unit 
vector parallel to the normal to M'^ at x, and ds is the usual volume d-form 
on M'^, induced from the Riemannian metric on W^^^ by the embedding 
M'^ C W^^^. Of course, we need to assume that the normal to M'^ is defined 
almost everywhere, i.e. that M"^ is smooth enough. 

Suppose additionally that the hypersurface M'^ is bounded, closed and is 
thus a boundary of a compact region N'^'^^ C M°'+^; we denote the volume 
jiV^+^l of iV'^+^ by vol {M'^) , and will call it the volume inside M'^. We will 
denote by D the collection of all bounded, closed hypersurfaces M'^, embedded 
in R''"'"^. We will not distinguish two elements of D if they are different only 
by translation. We add to D an element -k, corresponding to a manifold 
degenerating to a single point in M'^"'"^. 

Let h > he a real number. Consider the droplet energy functional 
$ on D, defined by 



The physical meanings of these functionals are the following: we think 
about the elements of D as droplet shapes, and the functional Wr represents 
the surface energy of a droplet. In the regime of droplet condensation the 
energy $ of a droplet consists of two parts, and in addition to the surface 
term we have also the volume term, which encourages the growth of the 
droplet. 

This functional has a local minimum at the point M^j^^ = -k, when the 
droplet M'^ degenerates into a single point, with = 0. For small droplets 
M'^ it is positive, and it takes negative values for very large Ms, when the 
volume term dominates over the surface term. Therefore $ has to have sad- 
dle points, and these saddle points are of importance for us if we want to 
understand the process of droplet condensation. The following theorem de- 
scribes the structure of the critical points of $. It turns out to be surprisingly 
simple! 

Theorem 1 Let the functional $ be given by with h > and the function 
T satisfying the properties given above. Then $ has precisely two critical 
points in D. One is a local minimum at M^j^ = -k. The other one is a saddle 
point 



$ (M'^) = Wr {M'^) - h vol {W^) . 



(2) 



sdl — 




(3) 
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Here f is a dilatation factor, while Wr is the Wulff shape, defined by 

Wr = dK<; with K< = e : Vn (x, n) < r (n)} . (4) 
The value $ {M^^^) is given by 

Note that the Wulff body is a symmetric convex bounded region in 
M'^^^. If at some boundary point the region has a unique support plane, 
then this plane is of the form 

(n) = {x G M^+^ : (x, n) = r (n) } 

for the corresponding n. Therefore 

VOl{Wr) = J^Wr{Wr). (6) 

With this remark the Theorem |I| follows easily from the Theorem ^ of the 
next subsection. 



2.2 Wulff minimizing problem. 

Let now Dg be the collection of all closed hypersurfaces M'^, embedded in 
M*^^^, and such that the volume vol (M*^) inside Af^ equals q. The Wulff 
problem consists in finding the lower bound of Wr over Dg : 

iVi ^JJq 

as well as the minimizing surface (s) Wt'^\ such that Wr (Wr''^^ = wi^\ if it 
exists. 

Theorem 2 The variational problem Q) has a unique (up to translations) 
solution 



[WriWr)) 

which is a scaled version of the Wulff shape, see 
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The paper |[1'2| ] contains a simple proof that Wr ^VTr'^^ j < Wr (M) for 

every M G -Dg. It proceeds in four steps, as follows. Let M = dN, = dK^. 

1) We observe first, that 

\N _L ^x<\ - \N\ 
WriM)> hm ' ^ ^ ' ^ . (8) 

Here for A,B C M'^+^ we denote by A+B = [z e W^^^ : z = x + y,x e A,y e B] . 
To see that (|]) holds, we note first that by continuity it is enough to check 
dl) in case when is a polyhedron. Then (§) follows from the remark 
that for a hyperplane L (n) , orthogonal to n, the sum L (n) + is a strip 
around L (n) with the width < r (n) , while the contributions of the lower- 
dimensional faces of to + eK^\ are of the order of at most e"^. 

2) Next, we apply Brunn-Minkowski inequality, which states that 

\A + B\> + |5| 

It implies that 

\N + eK<\ > (|iV|i/('^+i) +£|K<|i/(^+i))(^+i) 

> lA^I + {d + l)e\Nf/^'^-^^^\K<\^/^'^+^\ 

3) Without loss of generality we can restrict ourselves to the case q = 
vol (Wr) ■ Then we further have 

' ^ ^ ' ^ >{d + l) vol (Wr) . 

6 

4) On the other hand, (Wr) = {d + I) vol (Wr) , see (§). ■ 



The uniqueness of the minimizing surface is proven in |[ri|| . The proof is 
much more involved. 



2.3 Two-dimensional case 

In what follows we fix g = vol (Wr) and omit it from our notation. The 
specifics of dimension 2 is that the minimizer Wr of the functional Wr is not 
only unique, but also is stable in the Hausdorf metric. To formulate that 
statement we first remind the reader that the Hausdorf distance pn between 
the two sets A,C is defined as 

{A, C) = min {inf [r : A C C + Br] ,mf [r : C C A + Br]} , 

where B^ is the ball of radius r. 
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Theorem 3 Let 'y be a closed curve in M^, with the area inside it equals to 
that of Wr- Suppose that for < £ < 1 

Then there exists a shift 'y of 'y by a vector x = x(7) in R^, such that 

PH {l.W^)<C^e, 

where C = C (r) > 0. 

Of course, Theorem |^ imphes uniqueness of the minimizing curve of the 
functional Wt-. The Theorem |^ fails in higher dimensions. 
For the proof, see |DKS[|, Sect. 2.4. 



2.4 Maximizing problem. 

In a dual problem we again have a function 77 of a unit vector, but this time it 
is defined only over the subset A'^ = S'^ H M'j^"^^ of them, lying in the positive 
octant. We suppose again that the function is continuous and nonnegative: 
rj (■) > 0. We assume additionally that 

7] (n) uniformly as n ^ (9 A'^. (9) 

Let now G C be an embedded hypersurface. We assume that for almost 
every x&G the normal vector iia, is defined, and moreover 

e A^ for a.e. xeG. (10) 

Then we can define the functional 

V,(G)= / r7(njds,. (11) 

JG 

In analogy with the Section 2.2 we introduce the families Z),, g > 0, of such 
surfaces G as follows: 
G G 5g iff 

i) G splits the octant in two parts, with the boundary 9M+"^"^ belonging 
to one of them, 

a) the (((i + l)-dimensional) volume of the body Q {G) , enclosed between 
and G, equals q. In what follows we denote the volume of Q (G) by 
vol (G) . 
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For example, let / (y) > be a function on W^, non-increasing in each of d 
variables, and G [f] C be its graph. Then 

vol(G[/])= / f{y)dy, (12) 



so if Jjgd / (y) dy = q, then G [/] is an element of Dg, provided the function 

/ is sufficiently smooth. 

Our problem now is to find the upper bound of over Di : 

Vr^ = sup Vri (G) , (13) 

as well as the maximizing surface(s) G such that (K?) = Vrj, if the 
maximizer does exist. Note that the last problem differs crucially from (|7|), 
since here we are looking for the supremum. In particular, this upper bound 
evidently diverges if taken over all surfaces, and not only over "monotone" 
one, in the sense of ([T0|), unlike in the problem (|^. 

It turns out that there exists a geometric construction, which provides 
a solution to the variational problem ([T3|), in the same way as the Wulff 



construction solves the problem (|^). It was found in ||S2 
Let 

K> = {x G : Vn G A'^ (x, n) > r/ (n)} , (14) 

G, = d{K>). (15) 

Because of (P), the surface Gjj is a graph of a function, (y) , y E M.'^, i.e. 
Gn = G [/^] . Note, that the integral vol (G^) (see (p!^)) satisfies 

vo.(G,)^M|) (16, 

(compare with (|^)), provided both sides of ([T6|) are finite. 

Theorem 4 Suppose the integral vol (G^) is converging. Then the functional 
Vrj has a unique maximizer, V^, over the set Di. It is given by the dilatation 
of the surface (|T^j; 

r d + 1 



V,(G, 

// the integral vol (G) diverges, then = oo. 
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2.5 Facets 



In this subsection we will discuss facets - the fiat pieces of the surfaces Wr 



and Vri- 



In general, different functions r can result in the same Wulff shape Wr- 
There is a way to make this correspondence one-to one, by requiring addi- 
tionally that the function r satisfies the pyramid inequality, introduced in 
PS|| . Let us formulate it for the two-dimensional case. Let A, B, C be three 



points in R . Consider three unit vectors nAB,i^BC,i^cA, orthogonal to the 
respective sides of the triangle ABC. The pyramid (or in this case triangle) 
inequality is the requirement that 

\AB\ T {uab) + \BC\r (n^c) > \CA\ r (iic-a) 

for every triangle ABC. The pyramid inequality is equivalent to the property 
that every plane 

L= (n;A) = {x G M'^+^ : (x, n) = A r (n)} 

is a support plane to the (convex) body (A) (see (^)). In the case when 
T is the true surface tension function of some statistical mechanical system, 
this property is expected always to hold. In case of the 2D Ising model 
the triangle inequality holds sharply, see |pKS|| . Sharp pyramid inequality 



implies smoothness of the corresponding Wulff shape. For some general class 
of models the pyramid inequality was shown to hold in ||M1V1-SK|] . In what 
follows we will suppose it to hold. The validity of the pyramid inequality is 
equivalent to the convexity of the afiine extension f of r to M'^'^^ defined by 
f (x) = ||x|| r (^n^^ • As a result, the function r is differentiable a.e., and it 

has directional derivatives at any point n G S''* along any tangent direction. 
We will denote by r' (j^n) the derivative of the function r taken at the point n 
along the unit tangent vector to S'^ at n. The set of all such unit tangent 
vectors will be denoted by S^~^. 

We say that the function r has a cusp at the point uq if for every u &S^~^ 
we have 

t' (u) ^ -r' i-u) . (18) 

(Note that at any point of smoothness of r we have r' (z/) = — r' (— z/) .) If 
that is the case, then the intersection of the support plane (no;A) with the 
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Wulff body (A) is a convex closed subset of this plane, with nonempty 
interior. Such flat piece of the boundary d {K^ (A)) is called a facet. The 
shape of the facet can then be found by an analog of the Wulff construction 
d^), applied to the function r' (u) on S^~^. This analog is especially simple 
if there is an additional symmetry of the function r at no - namely, that for 
every u G S^'^^ we have t' (z/) = r' (— z/) . Note that (|l^) then implies that 
t' (z/) > for every u G S^'^^. In this case the above mentioned analog is just 
the Wulff construction itself, see ||M-S|| . Namely, the facet shape is given by 



F (r, no) = {x G M'^ : Vz/ (x, z/) < r' (z/) } . (19) 

All said above is applicable to facets of the surfaces K). The only difference 
is that in addition to bounded facets they may have unbounded ones, lying 
in ortants C (9M^|.'*'^. The shape of the bounded facets is given by the 
straightforward analog of (|19]). The unbounded facet in, say, xy plane (in 
d + 1 = 3-dimensional case) appears if the function t] has nonzero derivatives 
T]' (u) at the point n^ = (0, 0, 1) along tangent directions u G [O, |] = to 
the spherical triangle = S*^ fl M'^. In such a case the facet shape is given 
by the analog of (|1^: 

F {ri, n,) = {x G : Vz/ (x, u) > r]' (z/) } . (20) 



3 Statistical mechanics 

The relevance of the Wulff construction for the rigorous statistical mechanics 
lies in the fact that it describes the asymptotic shape of the droplet of one 
phase on the background of another phase in the regime of phase coexistence. 
The flrst result of that kind was obtained in the book [ pKS|| . In what follows. 



we briefly summarize it, as well as more recent results. We will talk only 
about the Ising model, though more general models can be studies as well. 

The rest of this section is organized as follows. We flrst recall some basic 
facts about the Ising model. We then formulate the results about the shape 
of the droplets of the (+)-phase in the (— )-phase, in the two-dimensional 
case and in higher dimensions. We conclude by explaining the relevance of 
these results to the theory of met ast ability. 
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3.1 Basics about the Ising model 

Gibbs states: At each site in Z*^"*"^ there is a spin which can take values —1 and 
+1. The configurations will therefore be elements of the set {—1, +1}^''^^ = 
Q. Given a E fl, we write ct{x) for the spin at the site x G Z'^"'"^. Two 
configurations arc specially relevant, the one with all spins —1 and the one 
with all spins +1. We will use the simple notation — and + to denote them. 

Observables are just functions on Q. Local observables are those which 
depend only on the values of finitely many spins. 

We will consider the formal Hamiltonian 

Hnia) = E ^(^)^(^) - ^ E ^(^)' (21) 

x,y n.n. x 

where /i e is the external field and a E Q is a generic configuration. We 
define, for each set A CC Z'^^-'^ and each boundary condition ^ G f^, 

HA,^,h{cr)^~ E '^i^My)-\ E ^(^)^(?/) - |E^(^)- 

x,y n.n. x,y n.n. 

x,yeA xeA,y^A 

We introduce the subset C as Qa^^ = {a E fl : a (x) = ^ (x) for x ^ A} 
We further define the partition function 

ZA,^,T,h^ E cxp(-/3i/A,e,,(a)), (22) 

where P = 1/T. 

The Grand Canonical Gibbs measure in A with boundary condition ^ 
under external field h and at temperature T is defined on fl as 

//A,^,T,.(<T) - I Otherwise. ^^^^ 

For each value of T and h, the Gibbs measure iJ.K{i),±,T,h with (±)-boundary 
conditions in the square box A(/) of size / converges weakly, as / oo, 
to a probability measure that we will denote by fJ^±,T,h- li h ^ 0, then 
f^-,T,h = fJ^+,T,h, so it will be denoted simply by iiT,h- h — the same is true 
if the temperature is larger than or equal to a critical value Tc — Tc{d + 1), 
and is false for T < T^., in which case one says that there is phase coex- 
istence. The critical temperature depends on dimension; (1) — 0, while 
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Tc{d + 1) > for d > 1. The measure /i+,T,o = f^+,T is called the (+)-phase, 
and /i_,T ^ the (— )-phase. 

In case A is a cubic box, A(/), we can also talk about periodic boundary 
conditions in A, which correspond to wrapping A(Z) into a torus. Then 
we can similarly define the configuration spaces ^a{i) = the 
Hamiltonian H^i)^per^h{<^) for a G ^a{1), the partition function ZA{i),per,T,h 
and the corresponding Grand Canonical Gibbs measure fiA{i),per,T,h in A (/) - 
with periodic boundary condition under external field h and at temperature 
T. 

Next thing we need to introduce is the notion of the Canonical Gibbs 
measures. To do it we first fix a real number p, —1 < p < 1. Now, for every 
finite box A we consider the subset ^^^^ C ^^a,^, defined by the property: 

^el^l^^$^a(x) = [p|A|], 

where [■] stays for the integer part. The canonical partition function ^ ^ is 
defined by the analog of p^), where the summation range ^a,£, is restricted 
to ^^Ag. The Canonical Gibbs measure Pa^th i^ defined by the analog of 
p3|), where ZA,^,T,h is replaced by j,^, and Qa,^ - by f^A^. The measures 
fi'A{i)perTh introduced in the obvious way. One checks easily that the 
measures Pa^th^ ^^A(i)perTh ^^^^ ^'^ depend on the magnetic field h, 
so we will denote them by /iA(;),per,T' A^a.^.t- 

One can say that canonical Gibbs measures are obtained from the grand 
canonical one by conditioning according to the event J^xga'^ (•^) ~ Ip \M] ■ 
Another collection of states we will need is obtained from the grand canonical 
Gibbs distribution by conditioning over the event |^j;gA'^(^)| — [P l^l] ■ 
They will be denoted by /^f^r/i- Alternatively one can introduce the set 
Of^ by considering only the configurations a with | X^ieA I — [P 1^1] ' 
the corresponding partition function j, ^ and then again write the analog 
of (H). 

For the expected value corresponding to a Gibbs measure p'" we will 
use the notation (/);;; = J fdp-], where ... stands for arbitrary sub- and 
superscripts. The spontaneous magnetization at temperature T is defined as 
m*(T) = {a (0))+,T- It is known that m*{T) > if and only if T < T^. 

Contours. With every configuration a one can associate a collection 7 (a) 
of surfaces in M'^"'"^, which are called contours. They are defined in the fol- 
lowing way. Let us consider the set *B [a) of all lattice bonds x, y such that 
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(J (x) (T (y) = —1. The surfaces in the collection 7 (cr) are formed by the col- 
lection *P (a) of (i-dimensional plaquettes, which are dual to the bonds in 
^ (a) . Each connected component of ^ (cr) is called a contour. 

Surface tension: The direction dependent surface tension is defined in 
the following way. For each vector n e S'^, consider the configuration ^(n), 
to be used as a boundary condition: 



+1, if {x, n) > 0, 
-1, if {x,n) < 0, 



where (•, •) is the usual scalar product in M^+^. The surface tension along the 
plane perpendicular to n is given by 

TT n = hm - log 24 

i-*oo (jA {I, n) ^A(0,+,T,o 

where A{l,n) is the (d-dimensional) area of the intersection of the hyperplane 
{x :(x, n) = 0} with the box A(Z). It is known that for each T < Tc the 
surface tension tt(-) is a continuous strictly positive and finite function. 



3.2 "The big droplet has the Wulff shape". Dimension 
2. 

The measure /x^^^^ ^ on configurations can be viewed as measures on con- 
tours. The theorem below explains what are the typical properties of these 
random contours, which explanation validates the title of the present sub- 
section. 

In the theorem below we consider the values of the density p smaller than 
the average value (cr (0))A(i)^per,T ~ m*{T). That forces the extra amount of 
the (— )-phase in the (+)-phase. As the reader expects perhaps, this extra 
amount gathers into one big droplet, which has approximately the Wulff 
shape. 

The rest of this subsection deals exclusively with the case of dimension 2. 

Theorem 5 Let T < Tc (2) . Then for every density p with m*(T)/2 < |p| < 
m*{T) there exist a sequence Di of subsets of configurations, Di C ^^ly 
which are typical, that is 

and which have the following properties: 
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for every configuration a & Di the set 7 (a) of its contours has precisely 
one "long" contour, T (cr) G 7 (a) ; all other contours are not longer 
thanKlnl, K = K{(3,p), 

the area |Int F (cr)| inside T (a) satisfies 

I |Intr((T)| - A/2| < Kl^/^ (In/)", 

where X = ""^ct)^' ' ^ = ^ (p) ' 

there is a point x = x (a) E A(Z), such that the shift of T (a) by x (cr) 
brings the contour T (a) very close to the scaled Wulff curve, defined 
by the Ising model surface tension r (see ^^): 



PH{r{a)+x{a),J^^^^^^lWr] <Kl'/'{lnir. (25) 



The first version of tlie above tlieorem was obtained in [PKS|| . Tlie ex- 
tension to all subcritical temperatures is due to loffe and loffe&Sclionmann, 
PI , PI , Tlie value 2/3 of the exponent in (^) is an improvement of the 
original 3/4 result, and is due to Alexander 

The restriction m*(T)/2 < \p\ in the theorem is needed because without it the 
droplet may prefer to shape itself into a strip between two meridians rather than 
to take the WulfF shape, see |S1|. This strip shape is, however, disadvantageous 



once the volume of the droplet has to be small. 



3.3 "The big droplet has the WulfF shape". Dimension 
> 3. 

In higher dimensions the status of the analog of the Theorem ^ is the follow- 
ing: it might be valid for very low temperatures, and it very probably fails 
for intermediate temperatures below the critical temperature. 

What is known to hold is a weaker version of this theorem, where the 
distance between random droplet and the Wulff shape is measured not by 
Hausdorf distance, but is understood in sense. To state the corresponding 
theorem, we will associate with every configuration a from Qa{i) a real valued 
function (t) on the unit torus T'^"^^, and we then compare this function 
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with the indicator function Ixk<> where XK^ C T'^+^ is the Wulff body, 
properly scaled. 

The function M„ {t) is defined as follows. We denote by ii the natural 
embedding of the discrete torus A(/) into T'^"'"^, the image of ii being the grid 
with spacing y. For t G T'^"^^ we define bi (t) C T'^+^ to be the ball centered 

d+l / 1 



at t with radius and let Bi (t) C A(/) be its preimage under ii. Then 



x&Bi(t) 



(The choice of the scale is somewhat arbitrary here; the essential re- 

quirement is that the size of Bi (t) should grow with I, but stay infinitesimally 
small compared with l'^'^^.) 

We are going to discuss the typical configurations of the measure j, ^ 

for /i = 0. In order to have a droplet, which, moreover, does not wrap itself 
around the torus, we have to put restrictions on p. The following one: 



m*(T) (^1 - 2 (d + 1) ~j< \p\ < m*{T) (26) 

is enough. The straightforward calculation shows that under such restriction 
we have to expect to see a droplet XK^ with 



^ _ d+l/ ^+1 rn*-\p\ 

^'WriWr) 2m* • ^^^^ 

To make a precise statement we first introduce for every subset A C T'^^^ 
the indicator (t) = | t ^ A'^ ' ^very function v in we 

denote by U {v, 6) its 5-neighborhood in (T'^^^) . 

Theorem 6 Let the dimension d+l be at least two, and the temperature 
T < Tc{d + 1) satisfies some extra technical restrictions (see below). Then 
in the large volume and for densities p satisfying the function M„ (t) 
typically has the Wulff shape: for every 6 > 

lim/i<5),,,.,^J-^^<^(-)e U {U{hK<+t,5)uU{-I,j,<_,„S))\ = l. 
Here X is given by ^ 
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The shifts by all t-s of the Wulff shape XK^ appear in the statement 
since the location of the droplet can be arbitrary. The droplet F (a) can be 
either of (+)-phase or of (— )-phase, which is why we consider two indicators: 
^XK< ^''^d —IxK< ■ Note that if a point t is such that the ball Bi {t) stays away 
from the boundary of the droplet F (a) present in the configuration a, then 
the value M^- (t) should be expected to be ±m* (T) , depending on the phase 
inside/outside the droplet, which explains the factor — i;^- 

The range of the temperatures for which the theorem holds is a segment 
^0,Tc {d + 1)^ , except may be at most countable set of temperature values. It 

does contain all temperatures low enough. The value is the so called "slab 
percolation threshold", which is believed to coincide with T^. 



For a proof, see |CeP|| , or the review paper [PoIV 



3.4 Facets of the WulfF shape and of the random droplet. 



The following is known about the surface tension function r (see (^^) of a 
model of statistical mechanics. When the dimension {d + 1) of the model is 
two, then the functions r are expected to be smooth, so they have no cusps, 
and the corresponding Wulff curves do not have facets (=straight segments). 
This is known to be the case for the 2D Ising model (from the exact solution). 
The methods of the book |pKS|| allow one to prove that for a general class of 



the 2D models (so called Pirogov-Sinai models) the surface tension function r 
is in fact analytic. The cusps do exist in the low-temperature 3D Ising model 
surface tension, so there the Wulff shape indeed has corresponding facets and 
thus resembles a crystal. This was proven in |]BMF[, and extended later to 



other models in ||M-S|| . Moreover, it is expected that in the 3D Ising model 
the so called roughening transition happens at the temperature Tr < T^, 
at which temperature the surface tension starts to be smooth, the facets 
disappear, so the Wulff shape becomes round. This, however, remains a 
challenging open problem. 

The Theorem ^ above provides some information as to how closely the 
random droplet, forced into the system by a constraint (p6D, resembles the 



Wulff shape. However, this result is not strong enough to permit one to un- 
derstand whether the big droplet (which is unique in a typical configuration) 
itself has facets, in some sense. Using the notation of the Theorem we can 
formulate the following hypothesis: 

Conjecture 7 Let the temperature T be low enough. Then the following 
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event has ^^^^-^ j, ^-probability approaching 1 as I ^ oo : 
There exist six distinct 2D planes Li = (a) C T^, i = 1, 6, two for each 
coordinate direction, such that the intersections Li (IT (cr) are flat facets of 
r (cr) . Namely, for every i 

i) diam {Li nT{a))>Ci (T) diam {\K<) I (see g^)), with Ci (T) ^ ^/2j3 
as T — > 0; 

n) (^ndfg^ > C2 (T) , wUh (T) ^ 1/2 as T ^ 0. 

Hi) The asymptotic shape of the facets L^nr (a) is given by the corresponding 

WulfJ construction, described in Sect. 2. 5. 

As T — i> 0, the Wulff shape tends to a cube. This is where the value ■y/2/3 in 
i) and 1/2 in ii) come from. 



3.5 Metastability and critical droplets 

In this section we describe rigorous results concerning the problem of conver- 
gence to equilibrium. Intuitively one should think about the situation when 
one observes liquid water at negative temperature - i.e. the supercooled wa- 
ter - and one is interested to know how the process evolves and how long 
does it take for water to freeze. 

Our playground will be again the Ising model in two dimensions. Now 
we need to supply it with the time evolution. For this we will use what is 
known as Glauber dynamics. It is a Markov process on the state space VL, 
whose generator, L, acts on a generic local observable / as 

{Lf){a) = Y,c{x,^){f{a^)-f{a)\ 

where is the configuration obtained from a by flipping the spin at the 
site X to the opposite value, and c(x, cr) is the rate of the flip of the spin 
at the site x when the system is in the state o. In words, one can say that 
the dynamics proceeds as follows: at every site x the spin cr [x] is flipped 
randomly, independently of all other, with the rate c(x,cr), where a is the 
current configuration. Common examples are Metropolis Dynamics: 

Ch{x,a) = exp{-p{A,Hh{(r)y), 

or Heat Bath Dynamics: 

Chix, cr) = [1 + exp(/3A,i7;,(a))]-^ . 
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Here (a)+ = max{a,0}, and AxHh{a) = Hh{a^) — Hh{a). The spin flip 
system thus obtained will be denoted by (o"|^/i.t)t>0) where ^ is the initial 
configuration at time t = 0. If this initial configuration is selected at random 
according to a probability measure then the resulting process is denoted 
by {(y^T h t)t>o- It is known that the Gibbs measures are invariant with respect 
to the stochastic Ising models. Moreover, 

^T,h;t ~^ f^-,T,h, (yT,h;t ~^ f^+,T,h, aS t ^ OO. 

We will be interested in the case when h is positive, though small. Then 
there is only one invariant state, fi+,T,h, so the state fJ'-,T,h is equal to ^+,T,h, 
and (Jxh t ^ f^+,T,h, as t — * oo. (One should intuitively think about the state 
cr^^.j for t small as the supercooled but liquid water, and identify the state 
fi+,T,h with ice.) We want to control the convergence of the temporal state 
cr^f^.-t to the equilibrium, fJ,+,T,h, and to see, if possible, that during some quite 
long time the state cr^i^.f looks very similar to the (— )-phase /i-,T, while after 
some time threshold it changes suddenly and looks quite similar to the state 
fi+^T,h- It turns out that all the above features can indeed be established 
rigorously. 

If one starts to simulate the above dynamics on a computer, then the 
picture observed would be the following: one would see that droplets of the 
(+)-phase are created in the midst of minuses, which droplets are there for 
a while, and then disappear. That process goes on for a while, until a (+)- 
droplet big enough is born. This one then starts to grow and eventually fills 
up all the display. One guesses rightly that the relevant object here is the 
droplet energy functional $, introduced in (|^). We need a two-dimensional 
version of it, so the set D consists of all closed curves m in M^. A trivial 
adjustment is needed in the volume (=area) term: 

$ (m) = (m) - m* (T) /i vol (m) . (28) 

It is due to the fact that because of the way the magnetic field h enters in 
the Hamiltonian (^), the gain in the energy per unit volume obtained by 
replacing the (— )-phase by the (+)-phase is given by m*{T)h. We remind 
the reader that the functional $ has unique saddle point, which is the curve 
msdi = msdi (T, h) = j^^Wr, with $ (msdi) = 2X;;^W. ( W^) , see (|), (|). 
Let us define the critical time exponent Ac = Ac(r) by 
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Theorem 8 Suppose that T < T^., h > 0. Let v he either the {—)-phase fi^^T 
or S{u=-}- (In fact, any v "between" these two states would go.) Then the 
following happens. 

i) IfO < A < Ac, then for each n G {1, 2, ...} and for each local observable 

f, 

n-l 

E (/ (a^,.,=exp{A/M)) = E (/) + ^ (^") ' (30) 

i=o 

where 

h, if) = lim 

(We stress that in the last relation we are using the Gibbs fields corresponding 
to the negative values of the magnetic field.) In particular, 

E (^^,.,=exp{A/M (0)) = -m*(T) + 0{h). (31) 

a) If \ > Ac, then for any finite positive C there is a finite positive Ci 
such that for every local observable f , 



Ktoxp{A/M)) - (/)t,.| <Ci||/||exp|-^| 



(32) 



The relation (|30|) imphes that the family of non-equilibrium states (■)^ ^.^^ , h > 
0, defined for every local observable / by 

{f)T,h^ = E (/ {(^T,h;t=cxp{\/h})) ) 

is a C°°-continuation of the curve {{■)-^T,h, h < 0} of equilibrium states. This 
is true for every < A < Ac and every as above. The states {■)xh-\ 
the metastable states we are looking for. The relations (|30|), (pl|) should be 
interpreted in the sense that before the time exp {Xc/h} our temporal state 
is still "liquid", while (^) means that after the time exp {Xc/h} freezing 
happens. 

This theorem was obtained in PS|| . Let us explain the heuristics behind it. 
It has two ingredients. The first one is that the transition to the equilibrium 
is going via creation of droplets of the (+)-phase. The second one is that once 
such a droplet is created by a thermal fluctuation, with the size exceeding 
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a certain critical value, it does not die out, but grows further, with a speed 
V of the order of h. (This second belief can be expected to be correct only 
in dimension 2.) Let us see how these two hypothesis can give us the right 
answer. As we already said, to get to the equilibrium we have to overcome the 
energy barrier, given by the functional $. Subcritical droplets are constantly 
created by thermal fluctuations in the metastable phase, but they tend to 
shrink, as is dictated by the energy landscape $. On the other hand, once a 
supercritical droplet is created due to a larger fluctuation, it will grow and 
drive the system to the stable phase. Since the minimal height of the barrier 
is $ (iTisdi) , one predicts the rate of creation of a critical droplet with center 

at a given place to be exp < — [ . 



by a factor of 1/3. The reason for that is the following. Note that we are 
concerned with an infinite system, and we are observing it through a lo- 
cal function /, which depends on the spins in a finite set supp (/). For us 
the system will have relaxed to equilibrium once supp (/) is covered by a 
big droplet of the (+)-phase, which appeared spontaneously somewhere and 
then grew, as discussed above. We want to estimate how long we have to 
wait for the probability of such an event to be close to one. If we suppose 
that the radius of the supercritical droplet grows with a speed f , then we 
can see that the region in space-time, where a droplet which covers supp (/) 
at time t could have appeared, is, roughly speaking, a cone with vertex in 
supp (/) and which has as base the set of points which have time-coordinate 
and are at most at distance tv from supp (/). The volume of such a cone 
is of the order of {vtYt. The order of magnitude of the relaxation time, treu 
before which the region supp (/) is unlikely to have been covered by a large 
droplet and after which the region supp (/) is likely to have been covered by 
it can now be obtained by solving the equation 



Comparing with (^), we see that we miss the correct answer, exp 





This gives us almost what we want: 




(33) 
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If we would know that v ^ Ch, then it further gives us that 



trel 



exp <i -- 



1 $ (msdi) 1 



T 



(34) 



The hypothesis that v ^ Ch seems natural in our context, yet it remains 
an interesting open question. Happily, it is clearly seen from the relation 
(p3|), that to conclude p^ ) from it we need much less. What we are able 
to establish in fSSj], and what is enough for our purposes, is the following 



weaker estimate: for every e > there is a Ci > 0, such that 



.>C,exp{-i}. 



4 Combinatorics 

In this section we consider applications of results of Section 2 to combina- 
torics. We discuss the typical shape of a big Young diagram and that of a 
random skyscraper (=plane partition). We make a detour in subsection 4.4 
and discuss there the famous MacMahon formula and the g- analog of the 
hook formula. 



4.1 Typical Young diagrams 

A partition p of an integer is an array of non-negative integers rii > n2 > 
... > rik > such that Yl'i^i^i — ^- ^e specified by the sequence 

{vk} of integers, with = / iff exactly / elements of p equal k. It can also be 
described by the monotone function 

oo 

Its graph G [(pp] provides a graphical description of p and is called a Young 
diagram. The same name Young diagram will be used also for the set Y = 
1/2) : 4>p (2/1) > 0, < 1/2 < 0p (2/1)} C M^. The number of the unit cells 
in Y will be called the volume vol {Y) of Y. The set of all possible Young 
diagrams will be denoted by y. The subset of Young diagrams with A^ cells 
will be denoted by y^- 
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Similarly, a plane partition S of an integer is a two-dimensional array of 
non- negative integers Uij, such that for any i we have nn > > ... > fiik > 

for any j we have nij > n2j > ... > Ukj > while again YlTj=i^ij ~ 
N. One defines the corresponding function 0^ (yi, 2/2) , ?/i,2/2 > 0, in the 
obvious way. The function (ps [yi, 1/2) is monotone in each variable. Its graph 
G [(ps] will be called a skyscraper. Sometime we will call by a skyscraper 
the partition S itself. Again, N will be called the volume vol (S) of S. A 
support, supp (S) , is just the support of the function 0^. The set of all 
(finite) skyscrapers will be denoted by S. If B C = {(1/1,1/2) : l/i,l/2 > 0} 
is a subset, finite or infinite, then by C 5 we denote the subset of all 
skyscrapers with support in B. 

Higher dimensional analogs of the above objects are straightforward. For 
example, we can talk about spacial partitions (4D diagrams), etc. 

Many more objects of a similar type can be defined. For example, one 
can put restrictions on how the steps of the stair G can look: they can 
not be longer than 3 units, and their heights can be only 1,2 or 5, say. The 
same freedom is allowed in 3D, and above. 

Let us fix the number N, choose the dimension and the kind of diagrams 
we are interested in, and consider the corresponding set of all these dia- 
grams. There are finitely many of them, so we can put a uniform probability 
distribution on Vjy. (Here, again, variations are possible.) The question now 
is the following: how the typical diagram from the family looks like, 
when N 00? 



The first problem of that type was solved in the paper |VKer|, see also 



yij |V^ , pVZ|| . It was found there, that the typical 2D Young diagram 



under statistics described above, if scaled by the factor ^l/-\/ivj , tends to 
the curve Cvk, given by 

More precisely, for every e > the probability that the scaled Young diagram 
would be within distance e from the curve ( p6D goes to 1 as — 00. 
The heuristic way to obtain (pB]) (and similar results) is the following: 
i) Let A = (ai, 02) , -B = (&i, 62) be two points in Z^, with ai < 61, 02 > &2. 
We can easily see that the number # {A, B) of lattice staircases, starting from 
A, terminating at B, and allowed to go only to the right or down, is given by 
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(^'"'(b^-Sr'"^)- Therefore one concludes by using the Stirhng formula that 

lirn ^_ln#(A,B) = /iKB). (37) 

Here n^s is the unit vector, normal to the segment [A, B] , the limit A,B^ 
oo is taken in such a way that the vector n^s does not change, and for 
n = (ni,?i2) , a = ^^^^^ , the entropy density function 

h{n) = - (niln—^ — + na In ^^^^ — ). (38) 
V ni+n2 ni + 722 / 

a) One argues that the number of Young diagrams of the area N scaled 
by VN, "going along" the monotone curve y = c (x) > with integral one, 
is approximately given by 



-P^^/ " -=F-=? ^lH^(^)fdA. (39) 



Hi) Assuming that indeed the model under consideration exhibits under a 
proper scaling some typical behavior, described by a nice smooth non-random 
curve (or surface) C, one comes to the conclusion that the curve C should be 
such that the integral in (0), computed along C, is maximal compared with 
all other allowed curves. 

In general case one is not able to write down the corresponding entropy 
function precisely. The only information available generally is the existence 
of the limit of the type of (|37|), by a subadditivity argument. It should be 
stressed that even when the variational problem for the model is known, the 
main difficulty of the rigorous treatment of the model is the proof that indeed 
it does exhibit a nontrivial behavior after a proper scaling. 

The above program was realized in [|VT| , |V2|| , see also [pVZ]| , for the 2D 



case described above and for some other cases. In a class of more general 



2D problems was studied. The first 3D problem was successfully studied in 
CKeP||. The method of the last paper can also solve the skyscraper problem. 



see 



One implication of these results is that the asymptotic shapes of various 
random combinatorial objects introduced above are given by the construction 
described in subsection |2.4|of the present paper, for corresponding rj function. 
For example, in the case of Young diagrams we thus have 
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Corollary 9 In the notation of the Theorem ^ the curve Cvk, given by the 
formula ^3di ) coincides with the curve K?, which results from the application 
of our construction 1^1^, (17^, to the function rj (n) = h (n) given by 

the formulas 

Of course, this statement can also be easily checked directly. 

Professor Elliot Lieb had suggested to me during our discussion in May, 
2000, that using the relation ( ^91) one can recover some results of Hardy- 
Ramanujan-Rademacher about the behavior of the function |3^jv| , i-e. the 
number of partitions of the integer A^. Namely, in the main order in it has 
to be 



where h is given by (0), (|38|) , and the functional Vh is defined by (0). The 
integral Vh (Cvk) is in fact easy to compute, by using our Theorem ^. Indeed, 
using its notation, we know from it that 



Therefore 



Vh (Cvk) = V^VjG^. (41) 



On the other hand, a straightforward computation shows that the curve 
Cvk intersects the line {x = y} at the point ^ ^^"^ ^ , ^ j , while the curve 

Gh intersects the same line at the point (^h ,h = 

(In 2, In 2) . Hence, 

Cvk = —Gh. (42) 

TT 

Comparing (^) and (pOD , we find that ^2/V/i (Gh) = VQ/tc, so in view of 

m 



Vh (Cvk) = TT 
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to be compared with Hardy-Ramanujan-Rademacher result that 



N\ 



exp < 71 



2N 



see, for example, ||An 



4.2 Facets of skyscrapers 

In this section we will discuss the problem of the typical shape of a random 
skyscraper, as described in the previous section. As was explained above, 
this shape (3 is a (non-random) concave surface in M.^, given by the dual 
Wulff construction of section (See ||Ke|| and especially [|CKeP|| for more 



details about this and related asymptotic shapes.) We will explain now the 
proof of the following: 

Proposition 10 The surface & has three identical flat pieces (facets); they 
belong to three quadrants (K+) • , = 1, 2, 3, which together form the boundary 
dE?^. Moreover, the flat curves d {& H (M^).) C are just the Vershik- 
Kerov curves Cvk (P^A properly scaled. 

Indeed, let rj {■) be the corresponding entropy function, defined on = 
5*^ n M.^. Its definition is the following (compare with (pT])): 

For n G consider the plane Ln = {x G M^, (x, n) = 0} . Let Sn be the 
family of all surfaces Q in with the properties: 

i) every Q is made from the unit closed 2D plaquettes of the lattice Z^, 

ii) the orthogonal projection of Q onto is a homeomorphism. 
Denote by Qn (one of) the surface(s) from Sn which minimizes the distance 
dist (Q, Ln) . Introducing now the quantity # {Bpt, Qn) as the number of sur- 
faces in En, which coincide with outside the ball Br of radius R centered 
at the origin, we define r/ (■) by 

ln(#(i?^,gn)) 
ry(n) = hm -— . 

R^oo 7lK'^ 

We will argue now that the derivative rj' {u) of the function rj (■) at the point 
ez= (0, 0, 1) in the direction u G M^, tangential to 5*^ C M^, is nothing else 
but h{i') , see (|37|). Once this is established, the claim about the shape of 
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the facets of the surface (5 follows by applying our dual Wulff construction, 



see dg). 

To see that rl (z/) = h{y) , let G be the vector defined by two 
properties: rig = ae^ + hv for some a, 6 real; and {n.Q,&:^ = cos 9. In words, 
the vector ng is a result of going from in the direction u hj a. distance 6. 
Since t] (e^) = 0, we have to consider the limit lime_^o ^^^^ = r]' (u) . 

Let 6 be small. Then the surface Q^, viewed from the z axis, can be 
represented by the collection of periodic lattice staircases s_i, sq, si, ... , 
with the spacing between two consecutive s-s of the order of (sin5)~^. If 
a surface Q contributes to # {Br, Qn) , then it is also represented by the 
sequence of lattice staircases s'_i, Sq, s[, ... , with the following properties: 
i) only finitely many of s'^ are different from the corresponding sf, 
a) each coincides with its correspondent, Si, outside Br, 
Hi) different s'j-s might touch each other and have common bonds, but they 
cannot cross each other. 

The removal of the third restriction makes the stairs independent, and that 
immediately provides us with the upper bound 

i {v) <h{v). 

To obtain the opposite inequality, let us introduce the modified functions 
/le (■) , which are defined similarly to the function h{-) (see (|37D), with the 
extra restriction that only those stairs are allowed which deviate from the 
segment [A,B] by at most e^^ /2. Consider now a collection s'_i, s'q, s'^, ... 
of lattice staircases, which satisfy the properties i), ii) from above and the 
following property: 

Hi') each stair s'^ can deviate away from its correspondent, Si, by at most 
(sin^)"V2. 

Clearly, every such collection defines a surface, Q (..., Sq, s'^, ...) , which 
contributes to # {Br, Q^) . That immediately implies that for every > 

hsine {^) < v' {^) ■ 

Since /isme (^) y h{h') as 9 ^ (see, for example, [ [MM-SK[j ), that proves 
the equality rj' {u) = h (z/) 



26 



4.3 Independent variables representation and some gen- 
erating functions 

In this section we will describe one technical device which can be used to 
prove some of the above combinatorial statements. 

We will call a (plane) partition strict iff all the entries of the (two- 
dimensional) array are distinct. In what follows, we will denote by 112 (N) 
(vrl (A^)) the number of different (strict) partitions of the integer N, and by 

(N) (7r| (A^)) the number of (strict) plane partitions of N. 

It is not difficult to see that the generating function 02 (x) of the sequence 
{■K2 ii),i = 1, 2, A^, ...} is given by 

00 ^ 

0^(-)=nT^- (43) 

1=1 

Indeed, after expanding each factor of (l43|) , we obtain the following expres- 
sion: 

ri,r2,r3,... 

where the summation is taken over all infinite sequences {ri,r2,r^, ...} of 
nonnegative integers, which have only finitely many non-zero entries. Each 
sequence of r-s defines the Young diagram by the relation (^). In the same 
way, the generating function of the sequence tt'^ (A^) is given by 

00 

02 {X) = l[{l + X'). 
1=1 

From the point of view of statistical mechanics the expression ( ^31) is 
nothing else but the grand canonical partition function of a one-dimensional 
system of non-interacting spin particles (i, I = 1,2, with (i = 0, 1,2, ... . 
The particles Q are not identical; the statistical weight for the l-th particle 
to be in the state k is equal to (x') . The parameter x plays here the role 
of chemical potential. Let us associate with every particle (i a volume v (Q) , 
so when the particle is in the state k, its volume will be Ik. The statistical 
weight Wx (C) of a configuration C = Ci) C2,C3' ••• ^'^^ given by 

(C) = a;^(fi)+'KC2)+^.(C3)+..._ (44) 
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The fact that the particles are not interacting is manifested by the absence 
of cross-terms in the exponent in {^). The product (02 {x))^'^ (C) defines 
a probabihty distribution on ^-s. This is what is called grand canonical 
ensemble in the language of statistical mechanics. 

Note now that if x < 1, then the distribution (^4|) allows only configura- 
tions C with finite total volume vol {() = v (d) + v (C2) + v (Cs) + ••• • That 
means that only finitely many of the values (1,(2,(3^ ■■■ can be non-zero. Note 
also that from every such configuration ( we can build a Young diagram Y {() 
in the following way: we first put (i blocks 1 x 1 in a horizontal array, we then 
put (2 blocks 1x2 next to them (here 2 is the height of the block), followed 
by blocks 1x3, and so on, until we use up all non-zero entries of (. If we 
place these blocks from right to left, we get a genuine Young diagram, Y [Q . 
The (normalized) distribution (|^) defines therefore for every x a probability 
distribution Pa, on 3^ = {Y} . 

With every grand canonical distribution one can associate a canonical dis- 
tribution, which in our case is the conditional distribution IP^ (■ | vol (■) = A^) 
on 3^, obtained from P^. by restricting it to the subset = [Y C y : vol (Y) = 
and normalizing to probability distribution. The crucial observation now is 
the fact that the canonical measure (• | vol (■) = A^) does not depend on 
X and is moreover a uniform distribution on y^- This is evident from (^4]). 
Therefore the question about the behavior of the typical diagram, asked in 
section 4.1, is the question concerning the properties of a certain canonical 
ensemble. 

At this moment the use of the jargon of statistical mechanics pays back. 
The stat-mechanical folklore immediately recognize our question as standard, 
and in order to answer it one should use the standard recipe: the principle of 
equivalence of ensembles! This general principle in our case claims, roughly, 
that the properties of a typical Young diagram are the same in the ensem- 
bles Px (■ I vol (■) = A^) and in the ensembles P^. More precisely, the typical 
properties of the Young diagrams under P^. (■ | vol (■) = N) , asymptotically, 
as A^ — » 00, are the same as those of the Young diagrams under Pa;(Ar), where 
the value x (N) satisfies 

Ep^^^, (vol(F)) = Ar. 

In words, the value x [N) has to be such that the mean value of the random 
variable vol (Y) under Pa;(7v) is the same as the (non-random) value of vol (Y) 
under P^^ (■ |vol(-) = A^) , i.e. A^. The advantages of passing to Px(Ar) are 
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evident: the random variables Q are independent under Pj., for every x. Of 
course, there is price to be paid: one has to prove first that the equivalence of 
ensembles indeed holds. But this is a standard, though sometime hard, work. 



The same can be done in the 
one ends up with collection of 



This is exactly the approach used in 
case of the strict partitions, see again [pVZ 
non- interacting particles Q, which however are taking now only values and 
1. More general cases can be treated by allowing the particles Q to interact 
via some local interaction, see IB 



4.4 More generating functions and the MacMahon for- 
mulas. 

Inspired by what was said in the previous section, one would like to repeat 
the same program for the plane partitions. The generating function 03 (x) of 
the sequence [N) is again known. It is given by the famous MacMahon 
formula (|54D below. Before discussing it we will make some preparatory 
comments. 

Let S be some plane partition of N, or a skyscraper. Then we can define 
a partition 

p = Il{S) (45) 

of N (or a Young diagram) , by just listing all the entries of the two-dimensional 
array S in the non- increasing order. One would like to study the inverse map, 
which to every partition p of N corresponds the preimage (p) , in hope 
to obtain in this way the method of counting plane partitions. The difficulty 
here lies in the fact that the preimages (pi) , (^2) might contain dif- 
ferent number of plane partitions, even when pi,p2 are partitions of with 
the same numbers of summands. The reason is that the partitions Pi may 
contain repeating integers, and the smaller is their number the bigger is the 
preimage, being maximal for strict partitions. 

We still want to implement this idea of counting skyscrapers via counting 
Young diagrams. So we will introduce a subclass of plane partitions, which 
we will call pedestals. To define them we first have to fix a complete ordering 
(n,ji) y («2,j2) y ■■■ y (ikjk) y ■■■ on the set Zl = : i,j > 0}. 

The choice of the ordering is not important, once it is red out from some 
strict plane partition. More precisely, the ordering (ii,ji) y (^2, j2) y ■■■ y 
{'i'kijk) y ■■■ will be called allowed, if the following is satisfied: for every 
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integer k the array 




k-l 




iik>l 
otherwice 



is a plane partition (of a corresponding integer). In particular, the first site 
(ii, ji) in every allowed ordering has to be the site (1, 1) • We choose for the 
ordering >- any allowed ordering, and it will be fixed in what follows. 

Let now S be some skyscraper. Using the ordering we can define 
another one, in the following way: 



Clearly, different skyscrapers S' and S" can lead to the same orderings 
i^s') ~ i^S") ■ In such a case we will call them equivalent, 5" ~ S" . 

Wc denote by O the set of all orderings, defined by skyscrapers: O = 



A skyscraper *P will be called a pedestal, iff every skyscraper 5" ~ ^ is 
higher than That means that if 5" = {nij} , and ^ = {riij} , then n.^ > n.^ 
for every pair (i, j) . In this case we will write ^ = ^ i'^s) — ^ {S) , and we 
call ^ a pedestal of S. 

A more direct definition of the pedestal, which also proves its existence, 
is the following: let ^ be some allowed ordering, which coincides with >- at 
infinity; in other words, O. The pedestal ^ (:>=) = {Tiij} is defined as 
follows: 

i) for every two sites ^ G Z^, with the pair directly 

preceding the pair under ^, we have 



ii) the skyscraper ^ (^) has zero height at infinity. 

In words, the height of ^ (^) at the site equals to the number of times 
the ordering ;^ disagree with the ordering >- once we proceed from to 
infinity, following the ordering ^ . (By definition, this number is finite.) 

The set of all possible pedestals ^ will be denoted by V. Note that by con- 
struction, the mapping ^ -w:^«p defines a one-to-one correspondence between 




{ys,SeS}. 




ff i^',J')y{^",J"), 



V and O. 
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We are now in the position to define a bijection, 

bst-.S^V X y. (46) 

Namely, to a skyscraper S we correspond a pair, Y) , where = ^ (5*) 
is the pedestal of S, while 

Y = p{S-V{S)), 

see (|45|) . Here we treat S' and ^ (S*) as real- valued functions, and S — ^ (S) 
is just their difference. By definition, — (S) is again a skyscraper, so 
reading its entries in decreasing order produces a Young diagram Y. 

The inverse map, {bst)~^ : V x y ^ S, is described as follows. To every 
ordering O and every partition p = {rii > n2 > ... > Uk > ...} we first 
assign an intermediate skyscraper S {'^,p) - by putting p along . More 
precisely, we define S {'^,p) to be the plane partition, which equals to Uk at 
the site G which comes k-th according to the ordering . Now we 
define 

ibstyH%p) = '^ + siy^,p). 

The bijection ( ^6|) and its generalization (^) which follows are in fact very 
close to Stanley bijection, |St|, which explains our notation. 

Note that if bst (S) = Y) , then by construction vol (S) = vol + 
vol (Y) . Therefore the generating function (x) of skyscrapers is just a 
product: 

03 (x) = 92 (x) X 0p (x) . (47) 

Here g-p (x) = ^2^=0 ^kX^ generating function of the pedestals: the 

number Ofc equals to the number of pedestals with volume k. 

Note now, that everything which was said above, about the set S of 
skyscrapers built over Z^, can be literally repeated for the subfamily of 
skyscrapers satisfying the restriction that their supports belong to a fixed 
finite subset B C Z^. (The only thing one has to be careful with is the 
ordering >~: one has to assume that according to it each point in B precedes 
each point in Z^\i?. That puts a restriction on the shape of i? : it has 
necessarily to have the shape of a Young diagram.) Their generating function 
will be denoted by (x) , the corresponding set of all orderings will be 
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a subset of O, while the set of corresponding pedestals satisfies = 
V n . Finally, the bijection hst is generahzed to 

6|: X 3^1^!, (48) 

where 3^'^' is the set of all partitions with at most \B\ nonzero entries. If we 
denote the generating function of the set of partitions belonging to 3^'^' by 
02'^' (x) , then we have 



1 

1=1 



The relation (|47| ) becomes 

sf {x) = 0P' {x) X QvB {x) . (50) 

Note that for finite B the set of the pedestals is also finite, and therefore 
the function g-ps (x) is just a polynomial. 

The next generalization we need is to consider skyscrapers 4>s{yi,y2) 
which are not defined everywhere on Z^. That is, we fix a finite subset 
A C Z^, and we consider functions 4>s{yi,y2) on Z^\A satisfying all the 
properties of the skyscraper function. All the objects of the previous para- 
graph can be introduced in this new context as well. (Again, we want the 
ordering >- to be such that according to it each point in A precedes each 
point in Z^ \A. This forces A to be Young- diagram- shaped.) 

Finally, we put both constructions together, considering a pair of Young 
diagrams A C B C Z^ , and the corresponding set S"^^ of skyscrapers defined 
on Z^ \ A, which vanish outside B. The ordering >- should be such that each 
point in A precedes each point in B\A, while each point in B precedes each 
point in Z^ \ B. The generating function for S"^^ will be denoted by {x) ; 
we have also orderings O^^, pedestals and the bijection 



The relation (|47| ) is further generalized to 

\B\A\ 



03 " {x) = Qv^B (x) 



n jhf' 

/=1 
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again with the polynomial generating function Q-pAB (x). 

There is another remarkable expression for the function g^^ (x) , but only 
in case when the diagram i? is a rectangle, B = m x n. U A = ^, then it is 
the famous MacMahon formula: 

m n ^ 

i=l j=l 

which is also valid for B = 1?\ : 



03(-)=nnT3^- (54) 

%=\ j=i 



For general diagram A C B ([53| ) is generalized as follows: 



beB\A 

Here the product is taken over all |-B\A| unit cells of the diagram B\A, 
while h {b) is the hook length of the hook of the cell 6 : it is the number of 
cells of the diagram B\Am the column of b below b plus the number of cells 



in the row of b to the left of b plus 1. The relation ( ^51) can be found in ||HG |. 

Comparing the two expressions, and (|^, we see that the following 
holds: 

Proposition 11 The polynomial g-pAs (x) satisfies 

I ^ n!r(i-^') ,,,, 

In particular, for the rectangle we have 

YY^^ (i - x') 

The fact that in (^) and (|57| ) we have complete cancellations, and the 
rational functions in the r.h.s. are in fact polynomials looks quite surprising! 
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In the limit as x ^ 1 we get 

i\B\A\)\ 



Q-pAB (1) 



]lb&B\A h (b) ■ 

This is the famous hook formula for the dimension of the irreducible repre- 
sentation of the symmetric group with |-B\v4| elements, corresponding to the 
diagram B\A. The relation (^) is its g-analog. It was proven in ||Ker |, |Kr|, 



though it was known earlier. Our proof via using ( p2D and (|55| ) seems to be 
the simplest known. 

Let us now discuss the higher dimensional situation. Consider the func- 
tion g^^"-^'^ (x) , which is the generating function of the spatial partitions, 
i.e. 3D arrays with obvious decay restrictions. Long time ago MacMahon 
had conjectured, in analogy with (^9]), (p^), that 

m n k ^ 

i=i j=i 1=1 

Later he started to doubt it, and then it turned out that this formula is not 



correct, see [|An|| . However, the higher dimensional analogs of the bijections 
(|I6| ) (|5TD as well as the product representation (|50D are always valid. For 
example, we have for the parallelepipeds 2x2x2 and 2x2x3 that 



2X2X2 ^ ^^"-^"^"^-^ , (58) 

where the polynomials 0-p2x2x2 (x) and 3^2x2x3 (x) are: 

gp2x2x2 (x) = 1 + 2x^ + 2x^ + 3x^ + 3x^ + 5x^ + 4x^ + 8x^ + 
4x^ + 5x^° + 3x^^ + 3x^2 + 2x^^ + 2x^^ + x^^ 



gp2x2x3 (x) = 1 + 2x^ + 3x^ + 5x^ + 6x^ + 12x^ + 14x^ + 

+ 25x* + 29x^ + 41x^° + 46x" + GOx^^ + 68x^^ + 86x^^ + 
+ 96x^^ + llTx^^ + 123x^^ + 141x^^ + 137x^^ + 144x^° + 
+ 140x2^ + 144x^2 + 137x2^ + Ulx^^ + 123x^^ + 117x^^ + 
+ 96x2^ + 86x^^ + 68x2^ + 60x^° + 46x^^ + 41x^2 + 29x^^ - 
+ 25x^^ + 14x^^ + 12x^^ + 6x^^ + 5x^^ + 3x^^ + 2x^° + x^^ 
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as one can find with bare hands. (In fact, the above computations of the poly- 
nomials gp2x2x2 (x) and 0^2x2x3 (x) were done together with Oleg Ogievecky, 
using Mathematica.) 

One then checks that the polynomial nf=i (l — 2;') is not divisible by 
0-p2x2x2 (x) , nor the polynomial Y[]=i ~ by gp2x2x3 (x) . That shows 
that, in some sense, there is no higher-dimensional analogs of the MacMa- 
hon formula. However, by (|58|), (^) we can still write down the relevant 
generating functions, having for example that 

41x^ + 63x^ + 104x® + 152x'^ + 230x^° + 

+ 171x* + 276x^ + 450x^° + ... . 

4.5 Independent variables representation for the skyscrap- 
ers. 

Now we can, using MacMahon formula 

oo oo ^ 

S3 i^) = UUT^^^i 

i=l j=l 

try to obtain the independent variables representation for the skyscrapers. 
Following the approach of the section we introduce the ensemble of non- 
interacting particles, Qj = 0, 1, 2, indexed by i,j > 1. Again, the particles 
(ij are not identical; the statistical weight for the particle at the location 
(i, j) to be in the state k is equal to {x'^~^^~^)^ . The volume v (Qj) associated 
to the particle Qj when it is in the state k, is {i + j — 1) k. The weight of the 
configuration ( of the particles Qj is given by the same formula (PI). 

Pushing the analogy with the section further, one should associate 
with every configuration ( the collection of solid hooks: if the particle Qj 
is in the state k, we have to take k hooks of the volume {i + j — 1) . One 
such hook is made from [i + j — 1) unit cubes; one can visualize it as having 
its vertex cube at the location , and going down and to the left till it 
hits the coordinate axes. The temptation is now to build from every such 
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collection of hooks a skyscraper, so that different hook collections will be 
transformed into different skyscrapers. This, however, seems to be impossi- 
ble. What is possible is somewhat more complicated construction. Namely, 
before building a skyscraper from a given collection of hooks, one has to saw 
each hook in two pieces along its bisecting plane. In this way we obtain 
twice as many blocks as we had hooks, of half-integer volumes. It turns out 
that one can then rearrange those blocks spatially, using only translations, 
in such a way that the resulting shape will be a skyscraper, and the corre- 
spondence thus obtained is one to one. This rearrangement is a special case 
of Robinson-Schensted-Knuth correspondence, see |[Kn|| or [[BKnf . It is not 
at all straightforward, so we will not reproduce it here. 



5 Conclusion 

In this paper we have described the explicit geometric construction, which 
predicts the asymptotic shape of some stat. mechanical and combinatorial 
objects. It is worth mentioning that the method presented should work in 
combinatorics whenever the underlying probability measure has certain lo- 
cality property, namely that the distant portions of the combinatorial object 
under consideration are weakly dependent. This locality property is in fact 
the key feature behind the results obtained in all the papers cited. It also 
holds for the corresponding problems of statistical mechanics, like the validity 
of the Wulff construction, and is crucial there as well. 
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